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Beschreibung 

Rechnergesttitztes Verfahren zur parallelen Berechnung des Ar- 
beitspunktes elektrischer Schaltungen 

Die rechnergestUtzte Simulation elektrischer Schaltungen hat 
bei der Entwicklung sehr grofier Schaltungen, also Schaltungen 
mit einer sehr groflen Anzahl von Elementen, eine iinitier groBe- 
re Bedeutung erlangt. Insbesondere bei der Entwicklung von 
Computerchips mit einer Vielzahl, beispielsweise mehrere hun- 
derttausend Transistoren, erweist sich eine serielle Bearbei- 
tung fiir die Ermittlung der Schaltungsgrofien durch einen 
Rechner wegen des zu grofien Zeitauf wands als unbrauchbar. 

In der WO 98/24039 wird deshalb vorgeschlagen, eine groBe 
Schaltung zu partitionieren und die Partitionen jeweils von 
verschiedenen Rechnern bearbeiten zu lassen. 

tJblicherweise werden bei der Berechnung zuerst der Arbeits- 
20 punkt, das heiiit die Potenziale aller Knoten als Grundlage 
fiir weitere Analysen, wie beispielsweise Transienten- oder 
Wechselstromanalysen, ermittelt . 

Zur parallelen Berechnung wird in U, Wever, Q. Zheng et al . : 
"Domain Decomposition Methods for Circuit Simulation" 
(Proceedings of the 8^^ Workshop on Parallel and Distributed 
Simulation, PADS '94 Edinburgh, Scotland, UK, S. 183-186, Ju- 
li 1994) und in U. Wever, Q. Zheng: "Parallel Transient Ana- 
lysis for Circuit Simulation" (Proceedings of the 29^^ Annual 

30 Hawaii International Conference on System Sciences, S. 442- 
447, 1996) eine Implementierung des Newton-Verf ahrens vorge- 
schlagen. Nachteiligerweise kann hier, aufgrund schlechter 
Konvergenzeigenschaf ten, nur bei Vorliegen hinreichend guter 
Schatzungen des Arbeitspunktes Konvergenz erzielt werden. 

35 Derartige gute Schatzungen sind in der Regel jedoch bei gro- 
fien Schaltungen schwer oder nicht zu bewerkstelligen . 
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Der vorliegenden Erfindung liegt deshalb die Aufgabe zugrun- 
de, ein rechnergestutztes Verfahren zur parallelen Berechnung 
des Arbeitspunktes elektrischer Schaltungen zu schaffen, das 
eine einfache, sichere und schnelle Berechnung des Arbeits- 
punktes der Schaltung gewahrleistet . 

Diese Aufgabe wird erf indungsgemaB durch das Verfahren gemali 
Patentanspruch 1 gelost. 

Bei dem Verfahren wird zur parallelen Berechnung des Arbeits- 
punktes das Auf ladeverf ahren verwendet, welches sehr gute 
Konvergenzeigenschaf ten aufweist. 

Durch das Auf ladeverf ahren werden, wie beispielsweise in H. 
Spiro: "Simulation integrierter Schaltungen" {2. Auflage, R. 
Oldenbourg Verlag Munchen, Wien 1990) beschrieben, die in ei- 
ner elektrischen Schaltung vorhandenen dynamischen Elemente, 
wie Kapazitaten und/oder Induktivitaten, dazu verwendet, urn 
tiber eine Pseudo-Transientenanalyse den Arbeitspunkt zu be- 
rechnen. Die dynamischen Elemente werden dabei schrittweise, 
von einem Wert Null ausgehend, auf einen im Prinzip beliebig 
hohen Wert, insbesondere auf den Wert "1" hochgef ahren. 

In einer Ausf uhrungsf orm nach der Erfindung konnen zusatzlich 
an geeigneten Stellen, beispielsweise an oder zwischen mehr^J 
ren Knoten dynamische Elemente, wie beispielsweise Kapazita-^B 
ten Oder Induktivitaten, mit vorbestimmten bekannten Werten 
zusatzlich eingebaut werden. Hierdurch konnen Schwierigkeiten 
bei der Berechnung aufgrund einer zu geringen Anzahl von dy- 
namischen Elementen vermieden werden. 

Hierzu konnen die zusatzlich in die Schaltung eingebrachten 
dynamischen Elemente von einem Anfangswert, beispielsweise 
gegen Unendlich, far den die Arbeitspunktberechnung trivial 
ist, schrittweise auf Null verringert werden, so dass wieder- 
um die ursprungliche Schaltung simuliert wird. 
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In einer bevorzugten AusfUhrungsf orm der Erfindung wird we- 
nigstens ein Knotenpunkt, also ein Verbindungspunkt von we- 
nigstens zwei Strompf aden, der Schaltung mittels einer Kapa- 
zitat mit einem vorbestimmten Wert mit einem vorbestimmten 
Potenzial verbunden. 

In weiterer Ausgestaltung der Erfindung kann aber auch an je- 
den Knoten Uber alle Partitionen eine Kapazitat angeschlossen 
werden, deren zweiter Anschluss jeweils an einem vorbestimm- 
ten Potenzial, beispielsweise an Masse liegt. Diese Vorge- 
hensweise hat den Vorteil, dass der Arbeitspunkt, also die 
jeweiligen Potenziale der einzelnen Knotenpunkte, ftir einen 
gegen Unendlich gehenden Anfangswert ftir die Kapazitaten, 
aufgrund einer Losung der dann expliziten Gleichung, ftir die 
Schaltung trivial zu berechnen ist. Uber eine schrittweise 
Veranderung des Wertes fur die Kapazitaten kann dann durch 
geeignete Neuwahl des Wertes flir die Kapazitaten die Schal- 
tungssimulation verandert werden, bis sich eine Berechnung 
des Arbeitspunktes der Schaltung fur einen Wert der Kapazita- 
ten gegen Null oder nahezu gegen Null ergibt. 

Durch diese Implementierung des Auf ladeverf ahrens auch ftir 
parallele Berechnung einer elektrischen Schaltung kann vor- 
teilhaf terweise eine sehr grofie Schaltung mit einer Vielzahl 
von Transistoren selbst bei einer geringeren Zahl von dynami- 
schen Elementen auf schnelle und einfache Weise berechnet 
werden. Fur die jeweilige Neuf estlegung des Wertes flir die 
Kapazitat mit dem Ziel, diesen gegen Null gehen zu lassen, 
sind verschiedene Vorgehensweisen denkbar, wobei als Ent- 
scheidungskriterium der Schwierigkeitsgrad der Berechnung des 
Arbeitspunktes des jeweiligen vorhergehenden Schrittes in 
Frage kommt . 

Der Arbeitspunkt , der sich durch Losung einer nicht-linearen 
Gleichung flir jeweils einen bestiinmten Wert der Kapazitaten 
ergibt, kann beispielsweise iterativ mit dem Newton-Ver f ahren 
gelost werden. Fur die Wahl des nachsten Wertes flir "C" kon- 
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nen dann die Anzahl der zur Losung notwendigen Iterations- 
schritte verwendet warden, bis ein vorbestiininter Wert fur "C" 
unterschritten wird. 

5 Weitere vorteilhafte Ausgestaltungen der Erfindung ergeben 
sich aus den Unteranspriichen • 

Die Erfindung wird nachfolgend anhand eines in der Zeichnung 
dargestellten Ausf ilhrungsbeispiels erlauter. In der Zeichnung 
10 zeigt die einzige Figur ein Ablauf diagramm, in dem die ein- 
zelnen Verf ahrensschritte des Verfahrens dargestellt sind, 

Nach der Zeichnung liegt als Ausgangspunkt 1 eine Simulatio 
einer elektrischen Schaltung, beispielsweise in der Schal- 

15 tungsbeschreibungssprache SPICE vor. Diese wird in einem er- 
sten Verarbeitungsprozess 2, wie beispielsweise in WO 
98/24039 als "Clustering-Verf ahren" beschrieben, partitio- 
niert, so dass sich einzelne Partitionen beziehungsweise Tel- 
le der Schaltungen, die moglichst mit gleichem Schwierig- 

20 keitsgrad berechenbar sind, ergeben. 

Selbstverstandlich sind auch andere anwendbare Partitionie- 
rungsvarianten, wie beispielsweise die in N, Frohlich, B. 
Riess, U. Wever, Q. Zheng: "A new approach for parallel simu- 
25 lation of VLSI circuits on a transistor level" (IEEE Transa 
tion on Circuits and Systems - I: Fundamental Theory and 
Applications, Vol. 45, No. 6, June 1998, Seiten 601 bis 613) 
beschriebene "Ratio-Cut-Methode" oder selbst eine willktirli- 
che Aufteilung "per Hand" denkbar. 

30 

In einem weiteren Arbeitsschritt 3 wird nun jedem Knotenpunkt 
der Schaltung, also einem Verbindungspunkz von wenigstens 
zwei Leitern beziehungsweise Strompfaden eine geerdete Kapa- 
zitat hinzugefiigt. Selbstverstandlich ist es auch denkbar, 
35 eine Kapazitat, deren zweiter Anschluss an einem vorbestimm- 
ten Potenzial liegt, hinzuzufugen, wobei sich sowohl Poten- 
ziale als auch die Werte filr die Kapazitaten jedes Knoten un- 
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terscheiden konnen. Aus Effizienz-, Stabilitats- und berech- 
nungsokonomischen GrQnden liegt jede Kapazitat luit ihrem 
zweiten Anschluss an dem gleichen Potenzial, beispielsweise 
Masse, wobei zudem aus denselben Grtinden fur alle Kapazitaten 
5 in Schritt 4 ein identischer hoher Wert CO gewahlt wird- 

Fur diesen Wert CO wird nun in einem weiteren Schritt 6 fur 
jede Partition beziehungsweise Teilschaltung der Arbeitspunkt 
berechnet, wobei die erforderlichen Kopplungswerte, also die 
10 Werte ftir die Kopplungspunkte beziehungsweise Schnittstellen 
benachbarter Partitionen, fur die Berechnung der Arbeitspunk- 
te benachbarter Partitionen ausgetauscht und einberechnet 
werden. 

15 Hierbei kann aus Ef f izienzgrunden vorteilhaf terweise eine 

Partition, welche im Folgenden als "Master" bezeichnet wird, 
die Steuerung des Auf ladeverf ahrens ubernehmen. Der Master 
bestimmt dann den Anfangswert CO fiir die Kapazitaten, wobei 
es auch denkbar ist, den Anfangswert extern, beispielsweise 

20 vom Benutzer, vorzugeben. Dieser Wert CO wird dann an alle 
ubrigen Partitionen, im Folgenden als Slaves bezeichnet, 
iibergeben. Nachfolgend wird sowohl in der Master-Partition 
als auch in alien Slave-Partitionen der Arbeitspunkt berech- 
net, wobei selbstverstandlich der voile Quellenvektor, der 
die in der Schaltung vorhandenen Energiequellen darstellt, an 
der Schaltung anliegt. 

Ausgehend von einem zu losenden Gesamtproblem 

30 wobei m die Unbekannten des Masters, 
p die Anzahl der Partitionen, 
Si die Unbekannten der Partition i 

bezeichnet, fuhrt das Auf ladeverf ahren hierbei auf das Diffe- 
35 rent ialglei Chungs system. 
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wobei 



D=diag (C^ . . .C, 0, . . . , 0) unci 
t die Zeit bezeichnet. 

Dieses Dif f erenzialgleichungssystem kann beispielsweise mit 
dem impliziten Euler-Verf ahren gelost werden, was auf die 
nicht-lineare Gleichung 

h 




fuhrt. Dabei ist die Losung ziim Zeitpunkt fur k = 0, 1, 
usw. und h die Schrittweite t*"*^ - t^. Diese nicht-lineare 
Gleichung kann beispielsweise iterativ mit dem Newton- 
15 Verfahren 



v-l 



iz)(xr-x^)] 



fur n = 0, 1, usw* gelost werden, 




Bei der parallelen Berechnung werden dann vom Master nur die 
20 Unbekannten "m" berechnet, wobei die Unbekannten "si" fur den 
Master Festwerte darstellen. Damit eine Losung dieses berech- 
net werden kann, werden von den Slaves bzw. Slave-Part itionen 
in jedem Iterationsschritt des Masters die korrigierten Werte 
berechnet und dem Master mitgeteilt . Hierzu muss von den Sla- 
25 ves das nicht-lineare Gleichungssystem 
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gelost warden, wobei fi und Di die entsprechenden Gleichungen 
und Matrizen fur den Slave "i" bezeichnen, Diese nicht- 
lineare Gleichung kann ebenfalls mit dem Newton-Verf ahren be 
rechnet werden, wobei nicht notwendigerweise alle Iterations 
schritte bis zur Konvergenz durchgeftihrt werden mussen. 

Die Ergebnisse der Slaves werden daraufhin in das Gleichungs 
system des Masters eingesetzt, woraufhin der Master den Aus- 
druck m^ll berechnen kann. 

In einem nachsten Schritt 7 wird ein neuer Wert "Cneu" fur "C 
vom Master bestimmt, wobei auch die Slaves Vorschlage inachen 
konnen. Urn das Ziel, namlich einen Wert fur "C" gegen Null 
beziehungsweise kleiner als einen vorbestimmten Wert e und 
damit eine Berechnung der urspriinglichen Schaltung zu errei- 
chen^ wird die Wahl des neuen Wertes "Cneu" fur "C" in Abhan- 
gigkeit von der Schwierigkeit der Berechnung des vorhergehen- 
den Schrittes mit dem Wert "Cait" fur "C" getroffen. 

Hierbei sind verschiedene Vorgehensweisen denkbar, beispiels- 
weise die Wahl von "Cneu" in Abhangigkeit einer Analyse der 
Anzahl der Iterationsschritte, die der Master zur Berechnung 
des nicht-linearen Gleichungssystems fur den vorhergehenden 
Wert "Cait" fur "C" benotigt hat: 




y2 n < ni 

C^i, fur ni ^ n ^ n2 
r, .2 fur n2 > n 



C =< 

^ neu 
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Wobei "n" die Anzahl der Iterationen des Masters und "ni", 
"n2" vom Benutzer vorgegebene Parameter sind. Anstelle der 
Halbierung beziehungsweise Verdopplung von "Cait" sind selbst- 
verstandlich auch andere Strategien zur Verkleinerung bezie- 
hungsweise Vergrofierung von "Cait" moglich. Zusatzlich kann 
auch die Anzahl der Iterationsschritte, die die Slaves bzw. 
Slave-Partitionen zur Losung ihres nicht linaren Gleichungs- 
systems benotigt haben, berucksichtigt werden, beispielsweise 
durch die Wahl von 



f 

n = max 



V 1=1 »=i 




15 wobei ni, siavej die Anzahl der Iterationsschritte des Slaves j, 
wahrend der i-ten Iteration des Masters bezeichnet. Hierbei 
ist anzumerken, dass die Wahl eines Wertes ftir "C" nur einen 
Einfluss auf die Effizienz des Verfahrens, nicht jedoch auf 
den Arbeitspunkt selbst hat. 

20 

Nach Erreichen eines Wertes von fur "C" kleiner oder gleich 
einem vorgegebenen Wert e wird die Berechnung an der Verzwei- 
gung 5 abgebrochen, wobei in einem letzten Schritt der Wert 
fur C auf "0" gesetzt werden kann, Hierdurch steht als Ergeb- 
25 nis 8 der Arbeitspunkt der ursprunglichen Schaltung fest un( 
kann iiber Ausgabeeinheiten, wie beispielsweise Bildschirm, 
Drucker oder ahnlichem, ausgegeben und/oder in einem Speicher 
als Grundlage ftir weitere Analysen der Schaltung zwischenge- 
speichert werden. 

30 

Durch das erf indungsgemafie Verfahren kann ein Arbeitspunkt 
einer sehr groBe elektrische Schaltung vorteilhaf terweise von 
einer Vielzahl von Rechnern oder Prozessoren parallel berech- 
net werden, wobei der Nachteil bekannter paralleler Berech- 
35 nungsarten, namlich fehlende Konvergenz bei unglinstigen An- 
fangswerten, vermieden werden kann* Der Auf wand fur die Be- 
rechnung hinsichtlich Iterationsschritten und schrittweisem 
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Setzen eines Wertes fur "C" wird durch die parallele Berech- 
nung aufgrund der Verteilung auf mehrere Prozessoren in einem 
vertretbaren Rahmen gehalten. 

Der so berechnete Arbeitspunkt dient dann als Grundlage fur 
weitere Analysen, beispielseise der Wechselstromanalyse einer 
Schaltung. 
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Patentanspriiche 

1. Rechnergesttitztes Verfahren zur parallelen Berechnung des 
Arbeitspunktes elektrischer Schaltungen, 

- bei dem die Schaltung in einem ersten Schritt in mehre- 
re Partitionen partitioniert wird, 

dadurch gekennzeichnet/ dass 

- zur parallelen Berechnung der einzelnen Partitionen das 
Auf laderverf ahren angewendet wird. 

2. Rechnergesttitztes Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, dass an we-^ 
nigsten einem Knotenpunkt der Schaltung ein dynamischessH 
Element (C, L) vorgesehen wird. 

3. Rechnergesttitztes Verfahren nach Anspruch 2, 
dadurch gekennzeichnet, dass an je- 
dem Knotenpunkt der Schaltung ein dynamisches Element (C, 
L) vorgesehen wird* 

4. Rechnergesttitztes Verfahren nach Anspruch 2 oder 3, 
dadurch gekennzeichnet, dass jeder 
Knotenpunkt der Schaltung mittels jeweils einer Kapazitat 
mit jeweils einem vorbestimmten Wert mit jeweils einem 
Potenzial verbunden wird, so dass ein Arbeitspunkt der 
modif izierten Schaltung berechnet werden kann, 

5. RechnergestUtztes Verfahren nach Anspruch 4, 
dadurch gekennzeichnet, dass 
an jedem Knoten einer Partition eine Kapazitat mit dem- 
selben Wert (CO) vorgesehen wird. 

6. Rechnergesttitztes Verfahren nach Anspruch 4 oder 5, 
dadurch gekennzeichnet, dass 
jeder Knoten einer Partition mittels einer Kapazitat mit 
demselben Potenzial verbunden wird. 
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7. RechnergestUtztes Verfahren nach Anspruch 4, 
dadurch gekennzeichnet, dass 

an jedem Knoten aller Partition eine Kapazitat mit dem- 
selben Wert (CO) vorgesehen wird. 

8. RechnergestUtztes Verfahren nach Anspruch 4 oder 1, 
dadurch gekennzeichnet, dass 
jeder Knoten aller Partitionen mittels einer Kapazitat 
mit demselben Potenzial verbunden wird. 

9. RechnergestUtztes Verfahren nach einem der Anspruche 4 
bis 7, dadurch gekennzeichnet, dass 
das Potenzial auf Masse liegt. 

10. RechnergestUtztes Verfahren nach einem der Anspruche 4 
bis 8, dadurch gekennzeichnet, dass 

- bei. geeigneter schrittweiser Veranderung des Wertes (C) 
der Kapazitat jeweils der Arbeitspunkt der Schaltung 
berechnet wird und 

- dieser Schritt solange wiederholt wird bis die Werte 
der Kapazitaten nahezu Null sind. 

11. Computerprogrammprodukt, welches in einen Arbeitsspeicher 
einer Rechenanlage geladen werden kann, mit einem Softwa- 
recode fur das Durchfiihren des Verfahrens nach einem der 
vorhergehenden Anspruche, wenn das Computerprogrammpro- 
dukt auf einer Rechenanlage lauft. 



12. Datentrager mit einem Computerprogrammprodukt nach An- 
spruch 11. 
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Zusammenf assung 

Rechnergestutztes Verfahren zur parallelen Berechnung des Ar 
beitspunktes elektrischer Schaltungen 

5 

Die Erfindung betrifft ein rechnergestutztes Verfahren zur 
parallelen Berechnung des Arbeitspunktes elektrischer Schal- 
tungen, bei dem die Schaltung in einem ersten Schritt in meh 
rere Partitionen partitioniert wird, bei dem zur parallelen 
10 Berechnung der einzelnen Partitionen das Auf laderverf ahren 
angewendet wird. 



15 Hauptzeichnung ist die einzige Figur. 
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